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We derive for the first time the growth index of matter perturbations of the FLRW flat cosmolog¬ 
ical models in which the vacuum energy depends on redshift. A particularly well motivated model 
of this type is the so-called quantum field vacuum, in which apart from a leading constant term Aq 
there is also an dependence in the functional form of vacuum, namely A(ff) = Ao + Si'{H^ — Hq). 
Since \v\ 1 this form endows the vacuum energy of a mild dynamics which affects the evolution 

of the main cosmological observables at the background and perturbation levels. Specifically, at 
the perturbation level we find that the growth index of the running vacuum cosmological model is 
7 a thus it nicely extends analytically the result of the ACDM model, ja ~ 6/11. 
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1. INTRODUCTION 

Over the last two decades, studies in cosmology 
strongly indicates that we are living in a spatially flat uni¬ 
verse that contains 4% baryonic matter, ^ 26% dark 
matter and ^ 70% some sort of dark energy (hereafter 
DE) endowed with large negative pressure. The DE com¬ 
ponent plays a vital role in the cosmic history because it 
provides the necessary theoretical platform toward de- 
scribiM the accelerated expansion of the Universe (see 
Refs. [3,13 references therein). 

Although there is a general agreement concerning the 
main ingredients of the Universe, there are different pro¬ 
posals regarding the underlying physics which triggers 
the late cosmic acceleration. In fact the unknown na¬ 
ture of DE, which challenges the foundations of theoret¬ 
ical physics, has given rise to a plethora of cosmological 
models. The majority of such scenarios is based either 
on the existence of new fields in nature or in some mod¬ 
ification of the theory of gravity (for review see Ref. i) 
with the present accelerating epoch appearing as a sort 
of geometric effect. 

An alternative avenue that one can follow in order to 
explain the cosmic acceleration as well as to overcome, 
or at least to alleviate, the various cosmological puzzles 
is to consider a runnin g v acuum A(i7). The idea was 
proposed some time ago [j-Q and it might be the seed of 
future important developments of cosmology at a more 
fundamental level - see e.g. no for a review. This 
point of view is quite general and can be applied to the 
entire history of the Universe, including inflation and the 
“graceful exit” problem in the early Universe [HI- In this 
context, we do not need to introduce in the analysis new 
fields in the analysis nor modify the theory of standard 
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gravity [General Relativity (GR)]. In this cosmic ideol¬ 
ogy the DE equation of state parameter w = Pde/pde 
is by definition identical to -1, but the vacuum energy 
density is a function of time. Notice that there is an ex¬ 
tensive old literature on purely phenomenological models 
in which the cosmological term is a function of time [H- 
[13 , and a more recent series of works in which there is 
an attempt to connect the time evolution with funda¬ 
mental aspects of quantum field theory (QFT) in curved 
spacetime, which include the aforementioned references 
on running vacuum and others such as [isl - l^ . In this ex¬ 
tensive body of literature the time-evolving vacuum has 
been phenomenologically explored as a function of time 
in various possible ways, and in the more formal QFT 
approach it has been mainly investigated as a possible 
function of the Hubble parameter. The latter is the basis 
for what we will refer to as the quantum field vacuum 
model, in which A = A{H) and on which we shall mainly 
concentrate here. 

To test the above cosmological possibilities, it has been 
proposed that the so-called growth index, 7, of matter 
perturbations [ 23 ] can be used as an observational tool to 
discriminate between modified gravity models and scalar 
field DE models which obey general relativity. Nowa¬ 
days the accurate estimation of 7 is considered one of 
the most basic tasks in cosmological studies. Not sur¬ 
prisingly it has become traditional to study, for each pro¬ 
posed cosmological model, its background expansion as 
well as the growth index of matter perturbations, as in 
this way one may get an impression of the main cos¬ 
mological and astrophysical consequences of the model. 
For example, it has been found that for those DE models 
based on GR and characterized by a constant equation of 
state parameter, the asymptotic value of the growth in¬ 
dex is 7 ~ %Z-^5 P^[ ~ [^ • Obviously, for the concordance 


AGDM model (w = —I) we recover the nominal value, 
namely 7 Ri 6/II. As far as the modified gravity models 
are concerned the situation is as follows. In the case of 
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the braneworld gravity of [s^ we have 7 ~ 11/16 (see 
[13, IsH - ISSIl '). for some f{R) gravity models it has been 
found that 7 ~ 0.415 —0.21z for various parameter values 
(see BUS) and finally for the Finsler-Randers cosmo¬ 
logical model Basilakos & Stavrinos [s^l have shown that 
7 Ri9/14. 

In this work, we wish to investigate the growth index 
of the running vacuum model A{H) = Aq-I-S;/ — Hq). 
The denomination of running is related to the fact that 
it can be motivated within the context of QFT in curved 
spacetime and specifically from the point of view of the 
renormalization group approach [J[l], see [i3| andji] 
for recent reviews and references therein. This run¬ 
ning vacuum model generalizes the traditional ACDM 
model at the background level and can be put to the 
test. Currently the value of the dimensionless free pa¬ 
rameter 1 / is observationaUy allowed to be in the ball¬ 
park of \i'\ ~ 0(10“^) [isl - l^ [23 - 1^ . From its theoreti¬ 
cal interpretation [namely, as being the coefficient of the 
/3-function of the running A(i?)] it is a natural value, 
which in addition fits in with the existing theoretical es¬ 
timates Hi- 

To the best of our knowledge, we are unaware of any 
previous analysis of this kind applied to dynamical vac¬ 
uum models and for this reason we consider that it can 
be of theoretical interest, and maybe we can extract also 
some practical consequences. 

The structure of the article is as follows. Initially in 
section 2 , we briefly present the background cosmologi¬ 
cal equations. The basic theoretical elements of the lin¬ 
ear growth are discussed in section 3, while in section 
4 we provide the growth index analysis in the case of 
the running vacuum. In section 4 we compare different 
7 ( 0 ) parametrizations and, finally, in section 5 we provide 
some discussion and finish with our main conclusions. 


2. BACKGROUND EVOLUTION 

The physics of the dynamical vacuum model under 
consideration is based on the renormalization group 
(RG) in QFT according to the aforementioned references. 
Within this framework, the evolution of the vacuum in 
the current epoch is given by 

+ -Hi), ( 2 . 1 ) 

where Aq = A(iJo) = S^IaoHI and v is provided in the 
RG context as a “/3-function which determines the run¬ 
ning of the cosmological “constant” (GG) within QFT in 
curved spacetime [ 3 ] . The value of v is estimated through 
the upper bound \i>\ < I/(127r) ~ 2.6 x 10“^, which is 
approximate and is valid as an order of magnitude. It 
ensues from assuming that the effective masses of the 
heaviest particles involved in the loops for calculatin g th e 
/3-function is of the order of the Planck mass [1, fTsLTl^l , 
but in general it is not completely fixed. If there is a 
large multiplicity of heavy particles at a grand unified 
scale below the Planck mass, it could as well stay of 


order 10 Notice that the observational bound on v 
depends on the particular implementation of the model 
(EID, namely, on whether e.g., the vacuum exchanges 
energy with matter or not. In the simplest cases the 
natural theoretical estimate yields v = 10 “® — 10 “^ (see 
Sola’s Ref. [ 3 ] for details) and in these cases v is generally 
significantly smaller than the original upper bound. At 
the same time usin g co smolo g ical data it has been found 
that 1^1 = which is in agreement 

with the aforementioned theoretical expectations. Let us, 
however, point out that in models in which matter and 
DE are self-conserved the observational limits are weaker 
and they tolerate the order of magnitude v ^ 10 “^ from 
the original estimate, as shown in the recent work [33| . 

Dynamically speaking, since j^j <C 1 in all theoreti¬ 
cal implementations, it is easy to check that prior to the 
present epoch the low-energy behavior of the model tends 
to the usual AGDM model, but it is by no means identi¬ 
cal (for a recent review see Refs. l3, IE] and references 
therein). 

Gonsidering Eq. m, let us now focus on the deriva¬ 
tion of the Friedmann equations. Such a proceduce is 
perfectly allowed by the cosmological principle embed¬ 
ded in the FLRW metric. Namely, the A term may per¬ 
fectly evolve with the cosmic expansion, meaning that 
ultimately evolves with the cosmic time, t, but in general 
it depends on an intermediate variable, which in our case 
is the Hubble function, H = a/a, where a{t) is the scale 
factor and the overdot denotes a derivative with respect 
to t. The corresponding generalization of the Friedmann 
equations reads: 

S^Gptot = SttGp^ + A{H) = , (2.2) 

SirGPtot = SirGPm - A(iL) = -2H - 3H‘^ , (2.3) 


where the total energy density is ptot = Pm + PA (with 
PA = A/SttG the vacuum component of it) and Ptot = 
Pm + Pa is the total pressure. For the matter-dominated 
epoch, and of course also in our days, (Pm, Pa) = 
( 0 , —Pa), where the dynamical character of the vacuum 
does not alter the usual equation of state that it satisfies. 
This is an important point to remark. In fact, this obser¬ 
vation explains why the dynamics of the vacuum entails 
a corresponding modification of the local energy conser¬ 
vation for matter at fixed G. The outcome is that matter 
must exchange energy with the vacuum in order to fulfill 
the Bianchi identity, and this translates into the follow¬ 
ing generalized conservation law involving both matter 
and vacuum energy densities: 


Pm T 3iLpm — Pa ■ (^-4) 

This equation is actually not independent of (12.21) and 
(Ell), and therefore, using any two of them, it is easy to 
derive the equation of motion for the Hubble rate: 

H+^H^ = AttGpa = . (2.5) 

From (12.5p and the vacuum model equation (12.11) , we are 
able to determine the explicit form of P as a function of 
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time: 


H{t) = Ho 


/ Had - ,, 

\ - coth 

V l-v 

[3 1 

^Ho\/{VIao - v){l -v)t 


(2.6) 


where fl\o = 1 — f2mO and Hq is the Hubble constant. 
Utilizing H = a/a the cosmic time, t{a), follows: 


t{a) 


2 

3Ho a/(Hao v) 


(2.7) 


Inverting Eq. we easily determine the scale factor 
a = a{t). Therefore, inserting Eg. (12.71) into Eg. (12.61) we 
arrive at 

E^ia) = = Hao + , (2.8) 

where we have rescaled 

(2.9) 

1 — iz 1 — 1 / 


Notice that for ^ = 0 all the above formulas correctly 
reduce to the standard ones for the ACDM, and the 
rescaled parameters become the ordinary ones, Hio —t 
HiQ. Moreover, whether in rescaled form or not the cos¬ 
mological parameters obey the standard cosmic sum rule, 
namely Umo + ^^ao = 1 = ^mo + ^2 ao- 

Regarding the matter evolution, from Eqs. (113) and 
(12.31) we find H = —AnGpm, and combining the latter 
with Eqs. dm and (1^ we obtain a differential equation 
for the matter density: pm + 3Hpm = 3i/Hpm- Integrat¬ 
ing it (using Pm = aHdpm/da) we find 

Pmia)= pmoa-^^^-'K (2.10) 


Notice, that pmo is the matter density at the present 
time (a = 1), and therefore Umo = Pmo/Pco, where pco = 
3 Hq /SttG is the current critical density. As expected, we 
recover the standard matter conservation law pm ^ a~^ 
only for i/ = 0. However, thanks to iz 0 we can have a 
mild dynamical vacuum evolution; see Eq. dm- 

Defining flmia) = Pm(a)/Pc(a) it is easy to obtain, 
with the aid of Eos. (12.101) and (12.8L 

o 

g(„) - 

For convenience we also define 


(a) 


E^a) 


1 - 1 / ■ 


( 2 . 12 ) 


Differentiating Eo. (l2.12p and utilizing (12.81) we find that 


Subsequently, upon substituting Eg. (12.101) into Eg. (12.41) 
and integrating once more in the scale factor variable, we 
are led to the evolution of the vacuum energy density: 


PA(a) 


PAO + 


^ PmO 
1 - 


- 3 ( 1 - 12 ) _ 2 


(2.14) 


Once more for i/ = 0 the cosmological solutions of the 
running vacuum model under study boil down to the con¬ 
cordance ACDM cosmology, and in this case pA = Pad at 
all times. 

Finally, the observational viability of the current vac¬ 
uum model has been tested previously, and in the most 
recent analysis provided in Gomez-Valent et al. [2^ it 
is found that (D^,!/) = ^82 ± 0.012,0.0048 ± 0.0032) 
(see also Table 1 in Ref. [2^). For the rest of the paper 
we shall take this result as the basis for our estimates. 
Notice that, due to the rescaling [see the first equality of 
(12.91) ] we have Qmo = 0.283 ± 0.012. Recall that for the 
concordance A model {i/ = 0) fimo = 0.291 ±0.011 which 
is in agreement with the recent Planck 2015 results 0. 


3. LINEAR GROWTH 


In this section we concentrate on the basic linear equa¬ 
tion that governs the evolution of the matter perturba¬ 
tions. Following [ 1 ^ 123,123 we write the following equa¬ 
tion for the matter density contrast D = 5pm/Pm- 

b + {2H + Q)b - {^TTGpm - 2HQ - q) D = 0, (3.1) 

where Q{t) = —paIP m- For a formal proof of this equa¬ 
tion in relativistic cosmology, see Refs. [13, Hlj. It as¬ 
sumes that the DE perturbations are very small and that 
the divergence of the perturbed matter velocity is also 
negligible. Obviously, the running vacuum energy still 
affects the growth factor through the function Q{t), and 
therefore, it affects the background evolution of the mat¬ 
ter perturbations. 

In this context, we can rewrite the homogeneous form 
of Eg. (13.II) in terms of the scale factor (using d/dt = 
aH{a)d/da) as follows: 


a 


2 


d'^D 

da^ 


^3a ± a 


dh\H 

dlna 


aQ\ dD 
H ) da 



3 

2(1 ^')b-m 


a dQ 


H~d^ 

2Q 

a dQ 

H 

H da 


(3.2) 


where 


dlnH dlnE 
dlna dlna 



- v)nmia) 


(3.3) 


dUm 


-3(1 - v)Q.m(,a) 1 


^m(^) 


Q(a) _ PA (a) 
H{a) H{a)pm{a) 


(3.4) 


dlna 


(2.13) 






































and 


a dQ 
H(a) da 


y "" 

--z/(l - v)Q.m{a). 


(3.5) 


Notice, that in order to derive the above expressions we 
have utilized Eas. H2.8|) - (l2.14|) . The growing mode solu¬ 
tion of Ea. (l3.2|) is written as (for more details see Ref. 

0 ) 


D{a) = Cia^E{a) F 



1 3 J_ 
2 ’ 2 ’ ^ 


3 

2 ’ 





(3.6) 

where ^ = 1 — v, Ci is an integration constant to be 
constrained by an initial condition, and F is the hyper¬ 
geometric function [40j |. 

Of course for the concordance ACDM model = 0 or 
Q{t) = 0] Eg. (13.21) reduces to the standard perturbation 
equation a solution of which is (see Refs. BEl) 


hVlmoF{a) dx 

^ 2 io x^F^[x) 


(3.7) 


4. GROWTH INDEX 

For any type of dark energy, a useful parametriza- 
tion of the matter perturbations is based on the growth 
rate of clustering . In our framework, the natural 
parametrization is^ 


where flA{a) = 1 — flmia). At this point it is interesting 
to mention that there have been many theoretical spec¬ 
ulations concerning the functional form of the growth 
index and indeed various candidates have been proposed 
in the literature. In this work, we phenomenologically 
parametrize 7 ( 0 ) by the following general relation (4^ 

7 (a) =70 + 7i2/(a) ■ (4.4) 

In other words. Eg. (14.41) can be seen as a first-order Tay¬ 
lor expansion around some cosmological function y{a). 
The following options have been considered in the liter¬ 
ature: u}{a) = lnf 2 m(a) (hereafter Tq parametrization: 
0 ), a(z) (hereafter Ti parametrization: 0 ) and 2 ; 
(hereafter r 2 parametrization: 0). Below we briefly 
present various forms of 7 (a) for the running vacuum 
model in these various parametrizations. 


A. To parametrization 

In this parametrization we use y{a) = uj = lnf2m(a). 
Obviously, for z 3 > I, namely tlm{a) —1 (or w —>■ 0) the 
asymptotic value of the growth index becomes 700 ~ 7 o- 
Steigerwald et al. 0| proposed a general mathemati¬ 
cal treatment [see their Eqs. (5)-(I2)]^ which provides 
compact analytic formulas for the coefficients 70 and 71 . 
These authors start from the fact that for a large family of 
dark energy models (including those of modified gravity) 
the linear differential equation of the matter perturba¬ 
tions takes on the form 

D + 2vhHD - AnGyLpmD = 0 . (4.5) 


with fim(a) defined in Eq. (12.121) . The exponent 7 is the 
so-called growth index (see Refs. [28l - [^ Idll 0) and it 
plays an important role in cosmological studies as we 
discussed in the Introduction. Inserting the first equality 
of (14.1|) into Eg. (13.21) we derive after some calculations 


df 

dlna 


+ 


Q 

2 + 1 + 


dhiFI 

dlna 


f + f 


3(1 — v)^rn 2Q a dQ 

2 ir~ H~d^ ■ 


or 


df 


+ 


dlna 

3(1 - v)a 


2 + 2(1-^)^^ +Y 


f + f 


9!/(I-i/)~ 
— bv -\ --- 


(4.2) 


(4.3) 


^ In tlie running vacuum model 112.Il l one may check from 112.fil l 
that, at large redshifts z 3> 1, ~ 1 — t'. Therefore, if the 

growth rate of clustering is modeled as a power law, then it is 
more appropriate to use /(a) ~ ClmQF, because for z 3 > 1 we 
achieve the correct normalization Qrn{z) ~ 1 . 


Naturally, any modihcation to the Friedmann equation 
and to the theory of gravity is included in the quantities 
vh and y. For the nominal scalar field dark energy which 
adheres to general relativity we have vh = y = ^^ while 
for modihed gravity models we get vh = ^ and yt ^ 1. 
Notice that for the latter cases we have = 54^(a) 

by definition. In the following we will determine the pre¬ 
cise relation between the generic coefficients and 

our vacuum parameter u. 

To start with, we expect that if we allow interactions 
in the dark sector, then, in general, the cosmological 
quantities are different from unity. Furthermore, 

if the matter component evolves differently from the 
usual power law a~^ one may expect that the quantity 
^m{a) = Pm{a)/Pc{,o) is slightly different from Vlm{a) as 


^ The methodology of Steigerwald et al. |44| can be applied to the 
framework of 7 (a) = n\^'^ ■ However, for the purpose 

of our study we restrict our analysis to A = f. We would like to 
point that Eqs.(3) and (5) of Steigerwald et al. |44| have a typo. 
Indeed one has to replace there the quantity 1 + uh with 2i'h- 
Notice, however that this typo does not alter our results because 
in our case the coefficient uh is a constant [see Ea. ll4.12ll ] which 
implies that Nn. = 0, see Eg. 114. fill . 
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defined in Eq. (12.121) . For example, in our case, due to 
Eas. (l2.12[) and the first equality of (14.ip one can write 
Eg. (14.511 as follows: 


df 

dlna 


2vh 


d\TiH\ 


dlna 


f + r 


ijlQ.r, 


(4.6) 


where jl = /r(l — v). Similar to [4^ let us transform 
Ea. (l4.6p as 


dw , d'y. ,„ dlnTd 3 

-—(7+w^)+e ^+2vh+-^ - 

aina au) alna z 


<^(1-7) = 0 (4.7) 


where we have set w = Innm(a)- Within the mathemat¬ 
ical framework of [l^ one finds 


70 = 


3(Mo + Ml) - 2('Wi + Ni) 


and 


7i = 3 
-2 


2 2Xi 3Mo 

M 2 + 2MiBi(l — 2 / 1 ) -I- Moi? 2 (l — 2 / 1 , — 2 / 2 ) 


(4.8) 


2(2 -h 4Xi + 3Mo) 
-62(2/1, 2/2) + -^270 +?{2 + N2 


2(2 -h 4Xi + 3Mo) 

The following quantities have been defined: 


X„ = 


d"(duj/dlna) 


and 


Xn = 


d^l^H 


duj^‘ 


duj^‘ 


(jj—Q 


Mr,. = 


uj—0 


d^jl 

dtM 


uj—0 


Hr, = 


d" (dlnid/dlna) 


dw" 


(4.9) 


(4.10) 


(4.11) 


UJ—0 


with = 1. Also Bi 2 are the Bell polynomials 
of first and second kind, namely Bi{yi) = 2/1 and 
62 ( 2 / 1 , 2 / 2 ) = 2/1 + 2 / 2 - According to Steigerwald et al. 
[Till [see their Eq.(lO)] the pair ( 221 , 2 / 2 ) is equal to ( 70 , 0 ) 
implying Bi(l - 2 / 1 ) = 1 - 7 o, 62 (^ 1 ,^ 2 ) = 7 o and 
62(1 - 2 / 1 ,- 2 / 2 ) = (70 - 1 )^- 
For the current running vacuum model, the quantities 
vh and p, can be easily identified on comparing Eg. (14.21) 
with Eg. (14.61) . and we find 


l/ff = l + — = 1+ -!^ 

2H 2 


(4.12) 


and 

p = l — v — 


4Q 2a dQ 

-=-— = l — L' — - - 

SrimH da rim 


(4.13) 

where in the derivation of the second equalities we have 
used Eqs. (lOl) and (13.5p . Then, based on Eas. (l2.13L 
(13.31) . (14.121) and (I4.13L we obtain 


{Mq,Mi,M2,Ni,N2} = {l-2^/-3l/^4^/, -4i/,0,0} 


{Xi,'Wi} = {X 2 , 62 } = {3(1 


3(1-1/) 


}, 
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FIG. 1: In the upper panel we present the asymptotic value 
of the growth index as a function of v (solid line). Notice, 
that the dashed curve corresponds to 7 a ~ 6/11. In the lower 
panel we show the relative deviation [1 — 7Aff/7 a]%. 


and thus 


6 -I- 3j/( 1 — 3i/) 
ll-3z/(4-h3i/) 


(4.14) 


_ -4iy + 8iy{l - 70 ) + {1 - 2iy - 3 t/^)( 7 o - 1)^ 

~ 2[2-h 12(1-i/)-h 3(1-2 j/-3i/2)] 

27 o +6(1 - vho - 3(1 - v) 

2[2-h 12(1-i/)-h 3(1-2i/-3i/2)]' ^ 

If we take the aforementioned fitting value v = 0.0048 
from [ 2 ^, we find ( 70 , 71 ) = (0.5496, —0.009). For v = 0 
we recover the ACDM pair ( 70 , 71 ) = (^,—0.00729) as 
it shoul d 1^ . Lastly, since 1 / is of order of 0(10“^) 
[13, HH, l23|| it is safe to neglect high-order terms of v 
from Eq. (14.141) . In that case the asymptotic value of the 
growth index becomes 


7oo ~ 7o 


6 -I- 3i/ 

11 - 12i/ 


6 

11 


. 35 

1 -I- 1 / 

22 


(4.16) 


where 700 ~ [7(a)]a„=i. 

In the upper panel of Fig.l we show the asymptotic 
value of the growth index as a function of v (solid 
curve), where 12 parameter lies in the theoretical interval 
[—which implies that 'jah & [0.5235,0.5678]. In 
the lower panel of Fig.l we present the relative devia¬ 
tion [1 — 'jah / 7a] % of the growth index with respect to 
7A ~ 6/11. We observe that for negative values of 1 / 
the asymptotic value of the growth index becomes less 


2 








































6 


than 6/11 (the opposite holds for positive values). This 
deviation can reach ±5% when we attain the aforemen¬ 
tioned theoretical upper bound oiv = ±l/127r ~ ±0.026. 
These features can be easily understood from the approx¬ 
imate formula on the rhs of Eq. (14.161) . As advanced in 
Sec. 2, there are DE models based also on Eq. (EH) and 
generalizations there (cf. Ref, 1 in which the obser¬ 
vational limits lie at the border of the upper theoretical 
bound, so in general we can say that ^ 5% corrections 
to 7 are conceivable and could perhaps be within reach 
in the future. For the cases in which Eq. EH represents 
a dynamical vacuum model in interaction with matter, 
however, the corrections are smaller. Using the latest ob¬ 
served value of V = 0.0048 ± 0.0032 provided by Gomez- 
Valent et al. [2^ (see also Table 1 in Ref. [2^) we find 
7 Ajt = 0.5496 ±0.0028, to be compared with 7 a ~ 0.5454 
for the concordance model, hence a < 1% correction. 


B. ri ,2 parametrizations 


Here we generalize the original Polarski & Gannouji 
work. In particular, changing the variables in 
Eg. (14.31) from a(z) to redshift = —(1 ± z)~^^] and 
utilizing f{z) = we find 


m 

'1 


(1 ± z)"f'\n{(lm) ± 37(1 - 

HI ± 


3, 

-±-(l-^.)UA±y 


--- bv -\ --- 


(4.17) 


where prime denotes a derivative with respect to the red- 
shift. For those y{z) functions which satisfy the restric¬ 
tion y{z = 0) = 0 [or 7(0 = 0) = 70]^, we obtain the 
parameter 71 in terms of 70, and v. Specifically, if 
we substitute z — 0 and 7'(0) = 7i2/'(0) in Ea. (l4.17L 
then we have 


7 i 


^mo - 3(1 - J")( 7 o 


where 


1)Uao +1 ~ 1(1 ~ ^ 

j/'(0)lnUmO 

(4.18) 



<dv{l - 


1-70 




Clearly, in the case of the usual A cosmology (iz = 0, 
= 0) the above formula boils down to that of Polarski 
& Gannouji [i^ for y{z) = z {£2 parametrization). Fur¬ 
thermore, based on the ACDM cosmological model and 
for y{z) = 1 — a{z) = (Pi parametrization), we also 


® For ri,2 parametrizations 70 = 7(z = 0 ) is different with that 
of Steigerwald et al. l 44 l . namely 70 = [7(0)]^^ ^ 700 (see 

section 3 A). 
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FIG. 2: Upper Panel: The evolution of the growth index of 
matter perturbations. The curves are as follows. The results 
for A{H) in the Toq parametrization are given by the solid and 
dotted lines. Moreover, the dashed and the dot-dashed curves 
correspond to the ACDM in the Fo.i parametrization respec¬ 
tively. Lower Panel: Comparison of the observed and theo¬ 
retical evolution of t he g rowth rate f(z)(Js{z). The growth 
data can be found in [sJI (see their Table 1). For the ACDM 
we use (Ts = 0.829 while for the running vacuum model we 
have (78 = 0.758 (see [H. 1^1. 


confirm the literature results (see Refs, [i^ , [3 [i^ ). 

Now, due to the fact that the function y{z) = z goes to 
infinity at large redshifts for the rest of the paper we fo¬ 
cus on the Pi parametrization which implies y'(Q) = 1. 
In this case one can easily see that 700 — 70 ±71 as long as 
z » 1. Therefore, plugging jq = - 71 into Eq. (l4.18l) 

£and using at the same time 700 ~ 11-121/ derive 

the constants 707 in terms of ULmO and v. For exam¬ 
ple, in the case of = (0.283,0.0048), we find 

(70,71) = (0.5636,-0.140), while for the concordance A 
cosmological model we get (70,71) = (0.5565, —0.110) for 
(U„o,iz) = (0.291,0). 

5. DISCUSSION AND CONCLUSIONS 

In the upper panel of Fig. 2 we present the evolu¬ 
tion of the growth index for the running vacuum model 
A(i7) in the Pg parametrization (solid line), and in the 
Pi parametrization (dotted line). In the same figure 
the dashed and the dot-dashed curves correspond to the 
ACDM in the Pg and Pi parametrizations, respectively. 
The comparison indicates that the growth index of the 
A (77) and ACDM cosmological models is well approxi- 
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mated by the Fq and Fi parametrizations. Specifically, 
we find that the corresponding relative deviations are 

[1 7 7f;V7f;^] 0.7% and [1 - 7r^V7f;^] 7 0.6% 

which means that essentially both parametrizations are 
equivalent, namely they provide the same growth index 
results. Based on the previous analysis, it is interesting 
to mention that the differences that we find with respect 
to the ACDM growth index are near the edge of the 
present experimental limits. Indeed, in a recent anal¬ 
ysis of the clustering properties of luminous red galax¬ 
ies and the growth rate data provided by the various 
galaxy surveys, it has been found that 7 = 0.56 ± 0.05 
and Omo = 0.29 ±0.01 [HU. Obviously, the prediction of 
7 for all our possible cases lies within Icr of that range. 

In the lower panel of Fig. 2 we show the growth 
data (solid points; see Ref.j^) together with the pre¬ 
dicted f{z)as{z) — for the running vac¬ 

uum K{H) in the Fq (solid line) and Fi (dotted line) 
parametrizations, and for the ACDM (dashed line, Fq 
and dot-dashed line, Fi) We observe that the running 
vacuum model reproduces the growth data well, in a way 
that is compatible with the ACDM model. Notice that, 
in order to obtain the above results for the concordance 
A cosmology we utilize = 0.829, while for the running 
vacuum model we use erg = 0.750 (see Refs. [H, 

Also, L>(z) is the growth factor normalized to unity at 
the present time. 

As we have already said in the Introduction the deter¬ 
mination of the growth index is important in cosmologi¬ 
cal studies because it can be used as a tool toward testing 
the validity of general relativity on extragalactic scales. 
For those dark energy models that adhere to GR and 
characterized by a constant equation of state parameter 
it has been found that 7 ~ [2^-[31|, while for the 

A(i7)±GR case we obtained 7 a„ ~ 11-121/ - Obviously, 
the growth index reduces to 6/11 for the usual ACDM 
model (w = — 1 or 1/ = 0). 

In the case of extended theories of gravity the situation 
is as follows. Recently, for the holographic dark energy 
models it has been found that the asymptotic value of the 
growth index is 7 ~ 4/7 For the branewor ld g rav¬ 
ity of Ref.fs^l we have 7 ~ II/I6 (see also Refs.|3oi. Isl- 
l35l|l. for some f{R) gr avity models it has been found that 
7 ~ 0.415 —0.2Iz (see [^[^1 and lastly for the Finsler- 
Randers cosmological model Basilakos & Stavrinos 
have shown that 7 Ri 9/14. Based on the aforementioned 
results if the derived value of 7 (based on the next gen¬ 
eration of surveys, cf. Euclid) shows scale or time de¬ 
pendence or it is inconsistent with 6/11 then this will be 
a hint that the nature of DE reflects on the physics of 


gravity. 

To conclude, in this work we have analytically stud¬ 
ied the growth index of matter perturbations for the 
FLRW flat cosmological models in which the vacuum 
energy density is a function of the Hubble parameter, 
namely A(iJ) = Ar\ + — Hq). In previous compre¬ 
hensive studies we have utilized such a dy¬ 

namical vacuum model in order to investigate the back¬ 
ground expansion and have carefully compared the dif¬ 
ferences with respect to the concordance ACDM cosmo¬ 
logical model. We believe that the combination of the 
works of Refs.fl^ [21I, together with the current arti¬ 
cle provide a rather complete investigation of the obser¬ 
vational status of the RG running vacuum model both at 
the background and perturbation levels. 

Within this framework, we have calculated for the first 
time (to the best of our knowledge) the asymptotic value 
of the growth index, which is given by jah ~ ii-iL - 
Obviously, the obtained formula analytically extends in 
a very clear way that of the ACDM model, 7 a ~ 6/11. 
In our study we have applied the two most popular 
parametrizations for the evolution of the growth in¬ 
dex: 7(z) = 7o ± liy{z), with y{z) = InDm(^) and 
y{z) = z/(l ± z), we have solved the problem analyt¬ 
ically and we have thus provided for the first time the 
coefficients 70 and 71 in terms of 0^0 and v. The com¬ 
parison shows that the above "/(z) parametrizations are 
practically equivalent and the corresponding evolution of 
the growth rate of clustering matches quite well the re¬ 
cent growth data. Finally, we have estimated the nu¬ 
merical corrections that the running A(7F) model could 
produce on the growth index as compared to the concor¬ 
dance model and pointed out that they could reach, in 
some cases, the level of a few percent, hopefully accessible 
in the future. 
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